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Abstract 

We present here a detailed derivation of an explicit spin-dependent expression for the bending 
angle of light as it traverses in the equatorial plane of a spinning black hole. We show that the 
deflection produced in the presence of the black hole angular momentum explicitly depends on 
whether the motion of the light ray is in the direction, or opposite to the spin. Compared to 
the zero-spin Schwarzschild case, the bending angle is greater for direct orbits, and smaller for 
retrograde orbits, confirming our physical intuition about the loss of left-right symmetry from 
a lensing perspective. In addition, we show that for higher spins, the effect is more pronounced 
resulting in tighter winding of direct orbits with respect to the axis of rotation, and a higher degree 
of unwinding of retro orbits. A direct consequence of this effect is a shift in image positions in 
strong gravitational lensing. 
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I. INTRODUCTION 



The deflection of light as it traverses through curved spacetime remains one of the few 
available observational tools to quantitatively study the geometry surrounding a strong 
source of gravitation. The study of gravity's effect on light, or gravitational lensing, began 
with Einstein's prediction in 1913, and the subsequent confirmation by Eddington in 1919, 
of the bending angle of light in the relatively weak gravitational field near the sun. The need 
for detailed study of gravitational lensing in both the weak and strong deflection regimes has 
become more compelling as the precision of our observational tools have grown many-fold 
in the past few decades. 

Be lld in g an g le calculated Q S ** Sc— Id and Kerr geometries show that 
as we approach the depths of the gravitational potential, the bending angle exceeds 2ir 
indicating that multiple looping of a light ray around the center of attraction is possible 
(see for example, [5], page 678). In strong deflection gravitational lensing, images formed 
as a result of this are referred to as relativistic images. In order to study these images, we 
need analytical expressions for light deflection in exact form and sometimes as perturbative 
series expansions. In this paper, we present a detailed analysis of deflection of light in the 
equatorial plane of the spinning, or Kerr, black hole. Indeed, the program was started by 
Darwin[l|, continued by Boyer and Lindquist [6] , Chandrasekhar[7] and many others since. 
The key difference between our approach and earlier work is that we obtain an explicit 
expression for the bending angle for both direct and retrograde motion. Our final result for 
the bending angle is expressed explicitly in terms of the black hole mass m. = GM/c 2 and 
its spin parameter a = J/Mc, where J/M is the angular momentum per unit mass of the 
black hole. 

Our goal is to continue to work in the same perturbative framework that was started in 
for the Schwarzschild black hole and apply it to the case of deflection of rays confined 
to the Kerr equatorial plane. Bending angles, and therefore the position and magnification 
of images, depend crucially on whether the light ray is traversing in the same or opposite 
direction to the direction of rotation. In the Schwarzschild case [8[ , we expressed the bending 
angle in terms of the invariant impact parameter b, and in fact found it useful there to 
use the variable b' = 1 — b c /b, where b c = 3\^3m 9 was the Schwarzschild critical impact 
parameter. For the equatorial Kerr case, we introduce a new definition for b' that includes 
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the non-zero spin of the black hole. We are able to define precisely the approach towards 
the critical impact parameter which is referred to as the strong deflection limit (SDL). The 
weak deflection limit (WDL), on the other hand, in both the Schwarzschild and the Kerr 
case are easily defined in terms of the impact parameter as the limit b — > oo. The detailed 



SDL and WDL series are presented elsewhere 

In Section 2, we set up the basic framework with definitions of the variables and sign 
convention. We also discuss the meaning of a number of limiting values as we cross-check 
our results with the zero-spin (Schwarzschild) case. We obtain a formal exact expression 
for the bending angle in terms of elliptic integrals of the third kind. A plot of the bending 
angle as a function of the impact parameter from the critical value all the way to infinity is 
presented. We believe that this plot has not appeared in literature before. 

The exact bending angle is plotted numerically for both cases to show a remarkable 
difference between the deflection of light on the direct side and the retro side. On the retro 
side, as the light ray traverses "upstream" it suffers a smaller deflection angle; smaller even 
than the Schwarzschild bending angle. This unravelling of the light ray as it tries to loop 
around against the spin direction supports our physical intuition. On the direct side, the 
bending angle is greater than the Schwarzschild value, resulting in a tighter looping of the 
light ray. Moreover, we show that as we vary the spin parameter a from low to high values, 
relativistic images are shifted inward on the direct side, and outward on the retro side. We 
derive series expansions of the bending angle in the strong and weak deflection limits, to be 



presented elsewhere [9], and show that the analytical results can be used to predict image 
positions and magnifications. 

II. THE EXACT BENDING ANGLE FOR KERR GEOMETRY 

We begin with the Kerr line element expressed in the Boyer-Lindquist coordinates (with 
9 = 7r/2 for the equatorial plane): 

ds 2 = g tt dt 2 + g rr dr 2 + g H d<p 2 + 2g t(f> dtdcj) (1) 
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where 



9tt(r) 



2m. 
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r 2 + a 2 + 
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(2) 
(3) 

(4) 
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where t = cr and m. = GM/c 2 (gravitational radius) with r physical time and M the 
physical mass of the black hole. From the Euler-Lagrange equations for null geodesies, we 
have 
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Couched in the quantity J z that appears even at this early stage, is the sign of the angular 
momentum. This will decide whether the light ray is traversing along the direction of frame- 
dragging or opposite to it. We will use the following notation to keep track of the two types 
of orbits as we go through the analysis. We define 



sb 



(7) 



with s = Sign(J z /S) and b is the positive magnitude. The parameter s is positive for direct 
orbits and negative for retrograde orbits as shown in Figure [U 

With a = J/Mc it is convenient to introduce the following notation: 

a Jc 
a = — 
m. 



(8) 



GM 2 

limiting ourselves to cases where < a < 1, with a = being the Schwarzschild limit and 
a — 1 being extreme Kerr. Orbits with s — +1 and b s > will be referred to as direct or 
prograde orbits; and those with s = — 1 and b s < as retrograde orbits. 
Rewriting the above expression, we have 
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FIG. 1: Sign convention for orbits as viewed from above. The spin axis points into the page in this 
figure. 



Next, letting u = 1/r and using 



we have 



and 
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Combining the above expressions and simplifying gives us 



[1 — 2m. u + a u 
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V#7 [1 -2m.u(l -a/6 s )] 2 
where the quantity B(u), a cubic polynomial given by 



2m. I 1 - — ! ii> - I 1 - — )u 2 + ^~ 



(9) 
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has a maximum of two positive roots and at most one negative root. We consider the case 
of one real negative root U\ and two real distinct positive roots u 2 and u 3 . The three roots, 



given in terms of two intermediate constants P and Q that allow us to line up the roots in 
the order u\ < u 2 < u 3 are given by 

P-2m.-Q 

Ui = 13 

4m.r 

U2 = - (14) 
P-2m. + Q 

u s = -; (15) 

4m. r 



By comparing the coefficients in B{u) to those in the original polynomial in equation (TT2 
we first obtain the following relationship between P and {a, b, s, r } 



2 

a \ / _ a 



r °^7 = ro 7T^t' (16) 

where a and 6 are positive quantities. This gives us the following relation between the critical 
parameters: 
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Note that we need to set ro = r sc , 6 S = 6 SC and P = 3m. to obtain the critical values. Note 
that when a = 0, r sc — > r c = 3m. in both cases. The critical value for ro is less than 3m. for 
direct orbits and greater than 3m. for retro orbits. 

Comparing the other coefficients of the cubic polynomial B(u) we obtain the following 
additional expressions: 

Q 2 -(P- 2m.) 2 1 



8m.ro 1 b 2 (1 - a/bsf 



(18) 



Q 2 = (P-2m.)(P + 6m.). (19) 

The intermediate variables Q and P can be eliminated by combining the above relations 
to yield a simple cubic equation involving the impact parameter and the distance of closest 
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approach given by 



I a 

r 3_ 6 2V V/_ r + 2m.o 2 (l- £ | = 0. 
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which can be solved to yield 



r 



26 

7! 



b 2 



cos 



/ 



COS 



a 



( « 2 V /2 



\ 



(20) 



This expression for the distance of closest approach and the invariant impact parameter be- 
comes important in both the SDL and WDL series expansions [9] completely in terms of the 
invariant normalized quantity b' . Also, note that this expression reduces to its Schwarzschild 
version (see e.g., equation (9) in when we set the spin parameter to zero. 

So, in the strong deflection limit, we have P = 3m. and the following expressions involving 
the critical quantities: 



3m. 
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(21) 



and 



{b sc + af = 27m 2 ,(b sc -a). 



(22) 



Equations (I2T]) and (|22|) agree with Chandrasekhar (eqn 82, p329). 

Combining these two equations gives us the following relationship between the critical 
values 

b 2 sc = 3r 2 sc + a 2 (23) 

To solve the cubic equation (|2"2"|) we need to consider the direct and retrograde motion 
separately. For direct orbits (i.e., s = +1) we have 



{b +c + af - 27m 2 {b +c + a) + 54am 2 = 



(24) 



the solution of which is given by 

6 +c = —a + 6m. cos 
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m. 



(25) 
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FIG. 2: The critical impact parameter in units of m, as a function of the spin parameter a. The 
Schwarzschild critical impact parameter is shown for reference. 



For retrograde orbits with s = — 1 we have 



0_ c - a) 3 - 27ml(b_ c - a) - 54am^ = 0. 



Solving this cubic equation yields 

6_ c = —a — 6m. cos 
We can combine the two cases and write 

b sc = —a + s6m. cos 
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3 V m. 

Note that the critical impact parameter is a function of the black hole spin. In Figure [21 we 
have plotted the critical impact parameter as a function of a. From equation (|22|) . we have 
the following: 

(b sc - a) (b sc + a) 2 
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Using the above, we obtain the following expression for the critical value of the distance of 
closest approach: 



2m. < 1 + cos 



- cos 



-i 



-sa 



m. 



(29) 



which agrees with equation 87 on page 330 in [7j]. Figure [3] is a plot of r sc /m. as a function of 
spin. For < a < m., the critical value of r$ lies in the following ranges: m. < r sc < 3m. for 
direct orbits and 3m. > r sc > 4m. for retro orbits. Note that r sc stays positive throughout, 
given that this is the radial coordinate. 
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FIG. 3: The critical radius in units of m. plotted as a function of the spin parameter a. The 
Schwarzschild radius, the radius of the photon sphere and the Kerr inner and outer horizon radii 
are shown for reference. 



In the extreme case when a = M, we have the following results: 

direct: s = +1 r +c = m. b +c = 2m. 
retrograde: s — — 1 r_ c = 4m. 6_ c = —7m. 

in agreement with equation 89 on page 330 in Q. 

The sign of the impact parameter carries the information about whether the light propa- 
gation is in the same or opposite direction to the spin. This appears to be the most important 
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difference in our approach, one that yields results pertaining to retrograde motion in a nat- 
ural way. We believe that details involving retrograde motion are harder to extract when 
the sign of a is used to keep track of the change in the direction of the black hole spin. 

There are two other polynomials that appear in equation We will denote the roots 
of 1 — 2m.M + a 2 u 2 = as u±, given by 

m. ± y/m 2 . - a 2 

u ± = - 2 (30) 

Using this notation we can now write the ratio of the polynomials in terms of partial fractions 

as 

1 - 2m.n(l - Us) _ C + | C_ , 31 > 
1 — 2m. u + a 2 u 2 (u + — u) {u~ — u) 

Solving for C + and C_, we obtain 



'm. + a/iti 2 - a 2 )2m. (1 - u s ) - a 2 



C + = v ' v •- / ' v ^ — (32) 

2a 2 \ m 2 — a 2 



and 



c = a 2 - (m. - y/m 2 - a 2 )2m. (1 - u 8 ) 

2a 2 \/ m 2 — a 2 



Let us now consider a light ray that starts in the asymptotic region and approaches the 
black hole, with ro as the distance of closest approach. It then emerges and reaches an 
observer who is also in an asymptotic region. Although the change in the <fi coordinate of 
the light ray appears in the equations of motion, we will consider the deflection of the light 
ray from its original path a. The change in <p and the bending angle a are simply related 
by a difference of tt. The bending angle can be written down in a straightforward manner 
by integrating the equation of motion for and subtracting 7r from it. The integrals therein 
can be separated into two parts as shown below to make the lower limit equal to the smallest 
root u\. The integrals are then readily expressed in terms of elliptic integrals of the third 



kind (see 



10], for example) to obtain the exact expression for the bending angle. These 



steps are shown below. 
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where Tl(n±, k) and H(n±, ip, k) are the complete and the incomplete elliptic integrals of the 
third kind respectively. The argument k 2 is defined through the elliptic integral as usual in 
the range < k 2 < 1. Note that in some references the variable is referred to as k 2 and 
in others simply as k. The order in which the arguments appear in U(n,ip,k) also varies 
between different references and in Mathematica. 

Remark: In Mathematica, the built-in mathematical function for the incomplete elliptic 
integral of the third kind EllipticPi[n, 0, m] is defined by 

-1/2 



I [l - n sin 2 6] 1 [1 - m sin 2 0] 
Jo 



dd 



and the complete elliptic integral of the third kind is EllipticPifn, m] = EllipticPi[n, ir/2, m] 
The other variables in the above expression are defined as follows: 
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It can be easily shown that in the limiting case when a — > 0, we have VL + = 1, Q_ = 
and n + = to give 

V 

= - 7r + iJ^[K{k)-F(i J ,k)}, (35) 

where K(k) and F(ip, k) are the complete and incomplete elliptic integrals of the first kind 
respectively. In addition, in the limit when m. = (i.e., h —>■ 0) we recover zero deflection 
as expected. 

We will use the following convenient notation: 

h = — - uj s = — and u> = — - (36) 
r b s mi 

with uj s taking on the appropriate sign for direct and retrograde orbit. Note that in the 
limit {uj s ,ujq — > 0}, we recover the zero-spin Schwarzschild case, and in the limit h — > 0, we 
have the zero-deflection flat metric limit. Further, we define critical parameters analogous 



to the Schwarzschild case in 



h sc = 1 + Us and r sc = (37) 

1 — UJ S llxr. 



We also define the variable 



^1-^1-3^(^1 (38) 
h sc \r J \l + u s 1 



"m.\ / 1 - u 8 \ q^o 1 _ 3m. 



with 

i- 3 r , 

r J \1 + uj s J r 
We have introduced these different quantities for the Kerr case, keeping in mind that they 

should go over to those defined in the Schwarzschild case smoothly when a is set equal to 

Jl 

zero. So, as shown above, as a — > 0, h sc — > 1 and we recover the definition of h' in [8J. In 
both cases, h — > at critical, and h — > 1 as r$ approaches infinity. 

From a lensing perspective, we are interested in impact parameters beyond the critical 
value (SDL) extending all the way to infinity (WDL). We define the dimensionless quantity 
b' as 

V = 1 - % (39) 
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where the insertion of the quantity s guarantees that the b' stays between and I. Note 
that this definition goes over naturally in the Schwarzschild limit: 

sb sc a ^o , 3V3m, 
l ~— * 

Now, some of the intermediate variables can be eliminated to rewrite all quantities in 
terms of h, h sc , ojq and u s as follows: 
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All of the above variables are to be substituted into the bending angle expression. We note 
here that the quantities ro, h, h sc , and u s depend on b, while uq = a 2 /ml is independent of b. 
Any quantity that has an "s" in the subscript takes on a negative sign for retro orbits. The 
bending angle itself stays positive since the sign of in the equations of motion is determined 
by the incident ray in the asymptotic region. In other words, as the ray approaches critical 
on the retro side, the overall deflection is still towards the black hole even though the extent 
to which it is bent is smaller than in the static case. Another way of looking at this is to 
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change the direction of the rotation axis from "into the page" in Figure [T] to "out of page" 
and look at the ray on the right side of the figure for which s — — 1 now; the bending angle 
is still positive. 

We now have an explicit expression for the bending angle (1541) via equations (14T)1) - (14"4"|) 
in terms of b', a and m., the remaining independent quantities. In Figures H] and \5[ we 
have plotted the exact bending angle as a function of b' for a = 0.5 and a = 0.99. Recall 
once again that the critical impact parameter depends on the spin parameter as well as on 
whether the light ray is undergoing direct or retrograde motion. As can be seen clearly from 
the bending angle plots, the bending angle is greater than the Schwarzschild value for direct 
orbits, and smaller for retro orbits. In addition for the higher spin (d = 0.99) case the effect 
is much more pronounced. Plots for other values of d look similar to this, and in the limit 
a = 0, all three curves merge as expected. 



111,1 


i i 1 i i i 1 i i i 1 i i 
Schwarzschild 






Kerr direct (c?=0.5) 






Kerr retro (a =0.5) 




'- 






\ 






\ 






\ 

V 






\\ 


















, , , i , 






) 0.2 


0.4 0.6 0.8 






b' 





FIG. 4: Exact bending angle plotted as a function of the normalized impact parameter. The strong 
deEection regime is towards the left of this plot as b' goes to zero, and the weak deflection regime 
is as b' — » oo. The Schwarzschild bending angle is shown for reference. 



One of the first remarkable features is the suppression of the bending angle on the retro 
side compared to the Schwarzschild case. On the direct side, light rays are bent more 
because of they are being swirled "downstream" with the spin, while motion of light rays 
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FIG. 5: Exact bending angle plotted as a function of the normalized impact parameter. The strong 
deEection regime is towards the left of this plot as b' goes to zero, and the weak deflection regime 
is as b' — > oo. The Schwarzschild bending angle is shown for reference. 



"upstream" results in smaller bending angles. In both the direct and retro motion, as in 
the Schwarzschild case, the bending angle exceeds 2ir, resulting in multiple loops and the 
formation of relativistic images. As a direct consequence of the difference between the direct 
and the retro side, as seen by an observer, these images will be offset to different amounts 
on either side. 

A plot of the bending angle divided by tt as a function of a for b' = 0.5 is shown in 
Figure M to illustrate the dependence on the spin. The solid line in this plot is the zero-spin 
Schwarzschild case shown for reference. 



III. CONCLUSIONS 



We have presented a detailed calculation of the exact bending angle formula for a light 
ray in the equatorial plane of a Kerr black hole. The frame-dragging in the Kerr geometry 
has a different effect on the light ray depending on whether it is traversing with or against 
the spin direction. The extent to which it is deflected from its original path is much higher 
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FIG. 6: Bending angle plotted as a function of spin for b' = 0.5. The Schwarzschild case is shown 
for reference. 



for direct orbits than for retrograde orbits. As a result, relativistic images on the direct side 
shifts inward as a increases, and outward on the retrograde side 9j. Higher order relativistic 
images can be studied by setting & = 2mr and using the solution in the lens equation. In 
order to demonstrate this, however, we first need to determine the series expansions of the 
bending angle in both the strong and weak deflection limit. We present these two series 
expansions in a companion paper [9]. A direct application of this to the case when the 
source, lens and observer are perfectly lined up is also presented in {9]]. 

Analytical solutions and results derived from them are crucial to our understanding of 
lensing images in general, and the Kerr geometry in particular. Even with our results here 
being exact, we have only considered the equatorial Kerr case. Our goal is to work with the 
full off-plane null geodesies in a similar manner to derive analytical solutions for predicting 
the behavior of a light ray as it spirals its way around the Kerr black hole before emerging 
towards the observer. 
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